INFINITE-ENERGY 2D STATISTICAL SOLUTIONS TO 
THE EQUATIONS OF INCOMPRESSIBLE FLUIDS 



JAMES P. KELLIHER 

Abstract. We develop the concept of an infinite-energy statistical so- 
lution to the Navier-Stokes and Euler equations in the whole plane. We 
use a velocity formulation with enough generality to encompass initial 
velocities having bounded vorticity, which includes the important spe- 
cial case of vortex patch initial data. Our approach is to use well-studied 
properties of statistical solutions in a ball of radius R to construct, in the 
limit as R goes to infinity, an infinite-energy solution to the Navier-Stokes 
equations. We then construct an infinite-energy statistical solution to 
the Euler equations by making a vanishing viscosity argument. 
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1. Introduction 

We develop the concept of a statistical solution to the Navier-Stokes 
(SSNS) or Euler equations (SSE) in the plane for an important class of 
velocity fields having sufficient decay of the vorticity at infinity to recover 
uniquely the velocity field from the vorticity. In particular, this class of ve- 
locity fields includes the important case of a vortex patch: a velocity field 
whose initial vorticity is the characteristic function of a bounded domain. 

Our starting point is the velocity formulation of a SSNS on a bounded 
domain given by Foias in [5j. (A highly accessible account of the theory of 
SSNSs is given in [6j, to which we refer often.) We adapt this formulation 
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slightly, of necessity changing the energy equality and using the same class 
of test functions as for homogeneous solutions in the plane (Section [6|). Our 
definition of an infinite-energy SSE is the same with the viscosity set to 
zero. We construct our infinite-energy SSNS by showing that it is the limit, 
in a special sense, of a sequence of statistical solutions on balls of radius 
R as R ^ oo. At the core of our approach is the expanding domain limit 
for deterministic solutions to the Navier-Stokes equations established in [9], 
extended to handle infinite-energy solutions. We then construct our infinite- 
energy SSE by making a vanishing viscosity argument. 

In the deterministic setting, the infinite-energy solutions that we consider 
correspond to an initial velocity lying in the space Em of [2] and [3]. A 
vector V belongs to Em if it is divergence-free and can be written in the form 
V = a + v', where v' is in L'^(M?) and where o" is a smooth stationary solution 
to the Euler equations whose vorticity is radially symmetric and compactly 
supported. (See Section [2] for more details.) A unique global solution to the 
Navier-Stokes equations exists and remains in the space Em for all time as 
long as the forcing has finite energy. The same can be said of solutions to 
the Euler equations if one imposes restrictions on the initial vorticity; for 
instance, that it lie in n This encompasses the case of a classical 
vortex patch — initial vorticity that equals the characteristic function of a 
bounded domain. (See also Corollary 12.31 ) 

A key parameter of any vortex patch is the total mass of its vorticity. 



Only when m = will the velocity field have finite energy (lie in L^), which 
excludes the case of a classical vortex patch. In a sense, m measures how 
infinite the energy is. 

When working with statistical solutions one would like to allow m to take 
on different values, because classical vortex patches that are nearly identical 
will typically have different values of m. Thus, we need to consider the 
spaces Em for all values of m simultaneously. 

We give a velocity rather than a vorticity formulation of our statistical 
solutions for several reasons. First, a vorticity formulation would require 
imposing higher regularity on the initial vorticity than required for solu- 
tions to the Navier-Stokes equations: it would be technically quite difficult 
to assume anything weaker than the initial vorticity lying in L^, as in [1]. 
Second, it would be hard to obtain convergence of the vorticity in the van- 
ishing viscosity limit, even with higher regularity of the initial data, without 
knowing that the velocity decays at infinity, and this does not come from the 
Biot-Savart law. Third, a start in this direction has already been made in 
[3] for time-independent solutions to damped and driven Navier-Stokes and 
Euler equations in the vorticity formulation. 

Constantin and Ramos do not specifically address infinite-energy solu- 
tions in [4]; however, their definition of a such a solution requires no change 
at all to encompass infinite-energy solutions and neither does their proof of 
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the vanishing viscosity hmit. Their construction of a stationary statistical 
solution to the Navier-Stokes equations as a long-time average of a deter- 
ministic solution to the damped and driven Navier-Stokes equations does 
assume finite energy. Nothing deep need be done, however, to extend their 
construction to allow infinite energy solutions (and so allow initial vortex 
patch data): one need only assume infinite-energy forcing. 

This paper is organized as follows: In Section [21 we define the function 
spaces in which we will work. We characterize the projection operator we 
will use to construct the initial velocities in Qpt in Section [3l In Section HJ we 
define weak deterministic solutions to the Navier-Stokes and Euler equations 
and give the basic well-posedness and regularity results for such solutions. 
The deterministic expanding domain limit of ^ is established for infinite- 
energy solutions in Section [5l We give the definition of a statistical solution 
to the Navier-Stokes equations in velocity form, for finite as well as infinite 
energy, in Section [H and construct an infinite-energy statistical solution 
to the Navier-Stokes equations in Section [71 showing that it is unique. In 
Section [8] we construct an infinite-energy statistical solution to the Euler 
equations using a vanishing viscosity argument. 

2. Function spaces and the Biot-Savart Law 

Let 

r^/j = the disk of radius R centered at the origin, 

with r^oo = I^^i and define the classical function spaces of incompressible 
fluid mechanics, 

Hr = H{Qr) = {ue {L^{nR)f: divu = 0, u-n = on d^n] , 
Vr = y(Ojj) = {n G {H^{QR)f: divn = 0, u = on , 
H = Hoo = H(M.^), V = Voo = F(M2). 
We endow Hr with the L^-norm. For Vr, we use the if^-norm: 

ll^lly« = ll«llL2(n«) + l|Vn||^2(^^). (2.1) 
Note, in particular, that 

Had we used the Poincare inequality to replace Equation (j2.ip with the 
equivalent norm that includes only the second term, as is normally done for 
a bounded domain, it would have introduced a factor of R in the right-hand 
side of the first inequality, preventing us from having a consistent norm with 
which to compare solutions on ^}r for different values of R. 

Our deterministic infinite energy solutions will lie in the space Em of [3]. 
A vector v belongs to Em if it is divergence-free and can be written in the 
form V = a + v', where v' is in L^(M^) and where a is a stationary vector 
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field, meaning that a is of the form, 

^= ("72 _^ P9ip)dp, ^ P9{p)dpj (2.3) 

with g in C^(M). The subscript m is the integral over all space of the 
vorticity, 

uj{v) = div^ — d2V^. 

Em is an affine space; fixing an origin, a, in Em we can define a "norm" 
by ||o" + v'll^^ = ||^^'|li^2(Q)- Convergence in E^, is equivalent to convergence 
in the L^-norm to a vector in Em- 

We will find it convenient to fix a choice of origin for Em as follows. For 
El we choose cti of the form Equation (j2.3p with uj{ai) supported in the unit 
disk and with 

/ uj{ai) = / uj{ai) = 1. 

We can then use am = mai as an origin for Em- 

Let Vo-i be a given fixed stream function for cji. As in in [10], Vo-i is 
radially symmetric with 

^PM = C2 + ^log\x\ (2.4) 
for all |x| > 1. For |x| > 1, |fTi(x)| = 1/ |x| by Equation ()2.3p . and 

l|cillui/-o \o ^ = 27r / —rdr + 27r / —rdr 

= 27r \og{R/{R - 1)) + TT [(i? - 1)^2 _ ^-2j ^ as i? ^ oo. 

Equation ([231) also gives AcJi = AV-^V^ri = V-^A^'^i = on Q!{ . 

The spaces ff/j and Em, or and Smnif^(M^), are the appropriate ones 
for initial velocities for weak deterministic solutions to the Navier-Stokes 
equtions, but for the Euler equations more regularity is required to obtain 
well-posedness. Rather than being as general as possible, we will assume that 
for deterministic solutions the initial vorticity lies in L°° for solutions on Q/j 
and in L^" fl L°° , for some pq < 2 for solutions on ^ . Slightly unbounded 
vorticities could be handled, as in [9], with little complication. This gives not 
only existence but uniqueness of the solutions. (The uniqueness of solutions 
for bounded initial vorticity is due to Yudovich [T2], as is the uniqueness for 
unbounded vorticities [13j.) 

Thus, we fix po < 2, and define the spaces 

Ym = {u(^Era-- uj{u) G L^o H , 

with "norm" 

II^IIy™ — W^We^ + \\^{'^ ~ ^ra)\\LPQnL°° 
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and 

Y{Qj^) = {ue H{Qr) n H^{nn) ■■ uj{u) G L°°} 

with norm 

ll'"llY(Cfl) = + ll^('")llLPonL°°(njj) ' ^ < oo. 

Because LP°{Qr) C L'^^Qr), using only the L°°-norm of uj{u) in the Y{Qfi)- 
norm would give a simpler, equivalent norm. We avoid doing this, however, 
for the same reason we avoided the use of Poincare's inequality in defining 
the V/j-norm in Equation (j2.ip . 

For statistical solutions in the whole plane, we do not want to assume that 
the value of m is fixed, so we must deal with larger spaces. For statistical 
solutions to the Navier-Stokes equations we will use 

E = IJ and = [j EmCi H^{R^) 

and for statistical solutions to the Euler equations we will use 

Y = IJ Y^. (2.6) 

mG-R 

E, Y, and E^ are function spaces, being closed under addition, with the 
norms 

Ikm + v\\e = M + \Ml2 ' II^IIy = II'"IIe + lk(^)llLPonL- ' 
Ikm + ''^IIei = Ikm + + ||Vu||^2(]j2) . 

These norms induce metrics on their respective spaces. Because H and V 
are separable, so too are E and E^. The space Y, however, is not separable, 
because L°°(]R^) is not. 

There exists a unique decomposition of any u in E, Y, or E^ of the form 
u = <7m + V, m (zM, V H. Given such a u, we define 

m{u) = m, a{u) = cr„(„). (2.7) 

Definition 2.1. We say that the support of a measure fj, on the function 
space X is bounded in X if 

supp^ C {n G X: < M} for some M < oo. 

If y is a subspace of X, we say that the support of a measure /i is {X,Y)- 
bounded if 

supp^ {u (zY: \\u\\-^ < M} for some Af < oo. 

That is, the support of /U lies in the subspace, but only its norm in the full 
space is controlled. 

Lemma 2.2. [Biot-Savart law] Let p be in [1,2) and let q > 2p/(2 — p). For 

any vorticity uj in LP(M2) there exists a unique divergence-free vector field u 
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in LP(M2) + L9(m2) whose curl is uj, with u being given by the Biot-Savart 
law, 

u = K*uj. (2.8) 

Here, K is the Biot-Savart kernel, K{x) = (l/27r)x^/ 

Proof. See the proof of Proposition 3.1.1 p. 44-45 of [3j. □ 

Corollary 2.3. For any vorticity to in (1 there exists a unique 

divergence- free vector field u in L°°(R^) whose curl is uj, with u given by 
Equation \2. <g|) . If uo is also compactly supported then u lies in E^, where 
m = f^^uj. 

Proof. By Lemma [2. 21 applied with p = 1, K *uj is the unique vector field in 
n L°° whose vorticity is to. But also, 

\\K *uj\\^oo < Wixn^K) *uj\\loc + 11(1 -xni)K *uj\\^oo 

< ilxni-?^|lLi II'^IIl°° + 11(1 - Xs^J^IIl-- II^IIli 

< C ll'^llLinL°° • 

Here, XA is the characteristic function of A. This shows that, in fact, K*uj is 
in L°° and is the unique such vector field. The last statement in the corollary 
follows from Lemma 1.3.1 of [3]. □ 

3. Projection operators 

Let Pvr : — > be restriction to ^Iji followed by projection onto Vr. Pvr 
is well-defined because as Hilbert spaces Vr is a closed subspace of H^-^{Uji), 
the space of all divergence- free vector fields in {H^ (^Iji))'^ , endowed with the 
inner product, 

We can describe Py^ explicitly by characterizing V^, the orthogonal comple- 
ment of Vr in H\- By definition, w is in VW" if and only if {w, v) fji (q ^ = 
for all V in Vr. Treating A.w as a distribution, integrating by parts gives 
{w — ^w,v) = for all v in Vr, where Vr = Vr H 'D{Q.r). It follows from 
this that w is in if and only if 

Aif — w = Vp (3-1) 

for some p in L^(r2) (see, for instance. Proposition Ll.l of |llj) and, of 
course, divw = 0. 

Now let u lie in H^^^{Qr) and let u = Py^u. Then w = u — u lies in 
so from Equation (j3.ip . 

An — u = Au — n + Vp in ^}r, 

divu = Ap = inQR, (3.2) 

u = on d^R. 

Equality is to hold in a weak sense in Equation (j3.2p . Since u is in H^, 
however, / = An — u + H is in H~^, which is sufficient to conclude that 
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li is in y and p is in L^. (See, for instance, Remark 1.2.6 of [H].) Also, 
the solution to Equation (j3.2p is unique because otherwise it would follow 
that —1 is an eigenvalue of the Stokes operator, —VA, where V is the Leray 
projector. But the Stokes operator is positive-definite, so all its eigenvalues 
are positive. 

The estimates involving the operator Py^ are hard to prove directly using 
this characterization. It is simpler to employ an approximate projection 
operator \Jji, and use the fact that projection into Vr gives the closest 
vector field in to the vector field being projected. (The same idea is used 
for projection into Hji in 17].) 

To define we need two cutoff functions, and h^. 

Let (fi in C°°{Qi) take values in [0, 1] and be defined so that = 1 on 
^1/2 and so that (^i and Vfi are both zero on OQr. Let (Pr{-) = ipi{-/R). 
Observe that ipR and VipR both vanish on B^Ir. 

Let g in C°°([0,3/4]) taking values in [0,1] be defined so that g{0) = 
g'{0) = and g = 1 on [1/2,1]. Then define Hr in C°°(Oi?) by HrIx) = 
g{R — \x\) for points x in il/j \ Hr^i and Hr = 1 on 0,r^i. Observe that 

\\hR\\c, < Ck, (3.3) 

k = 0,1, . . . , for constants Ck independent of R in [1, oo). Also, Hr = and 
VhR = on di^R. 

Definition 3.1. Define \Jr:K^ ^ Vr by 

VRiu) = v^(/i^(v.^ - ^.^(ii))) + v^(v?^v.) 

for u = am + V in Em- Here, Vd is the stream function for v chosen so that 
Iur ^v = Oon d^R. 

Lemma 3.2. Py^ maps continuously onto Vr with 

\\u - PvruWjji^^^^ < C ||u - criu)\\Hi(^^^\^^^^) + \m{u)\ (3{R), (3.4) 
where 

= \Wi\\HHnR\na-i) ^ as ^ oo, 

and 

\\PVr{u - am)\\vji < \\U - CTmWv , (3.5) 

llPv-flWlly^ < II-uIIei + C \m{u)\ . (3.6) 

Proof. That Py^ maps onto Vr is clear, and it is continuous because the 
restriction and the projection operators are continuous. 

To prove Equation ([33]), let u = am + v in Em n H^iM."^). Then 

\\u - ^RuWRHnu) < Ikm - U/?(Jm||Hi(nH) + " ^R^\\m(nR) ■ (3-7) 
It follows from the proof of Lemma 4.2 of [9] that 
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Also, letting ip = ip^^^ — ipa,n{R) and using Equation p.3p . 
\\am - ^Rcrm\\m(nR) = ll^'^V' - ^^{hRip)\\minR) 

< 11(1 - hR)am\\H^nR) + 

< 111 - ^rIIci hni\\m(nj,\nj^_,) + l|V/iij||ci M\m(nR\nn-i) 

< C \Wm\\m(nR\nn_,) + C M\m(nR\nii_,) 

Because \ i^R-i has width 1 and vanishes on its outer boundary, we 
can apply Poincare's inequality with a constant that is independent of R to 
give 

IIV'llL2(nfl\Qfi_i) ^ ^ II WliL2(Qfl\nfl_i) = \Wm\\L2(nn\nR_-,) ■ 

Thus, 

\Wm - U/J(Tm,||/^i(n^) < C lkm||_f/i(nB\Qfl_i) = \m\ Ikl ||//i(njj\Qfl_i) > 

which vanishes as i? — > oo by Equation (j2.5p . This gives Equation (j3.4p for 
Ur. Projection into Vr gives the closest element in Vr, so Equation (j3.4p 
holds for Py^. 

Equation (|3.5p follows easily: 

||PVfl(-u - crm)\\Y^ < \\u - amWvR ^ 11^ ~ 0-mlly ' 

the first inequality holding simply because Pvr is an orthogonal projection 
operator. 

To prove Equation (igTB]) . let u = am + v in £;„^ n H'^(M?). Then 

||PyflU||^^ < llPy^a^lly^ + llPVfl^^lly^ < Ikmlly^ + Ibllv^ 
< Ikmlly + llt^llv = l"i| Ikilly + W^Wv = C\m\ + \\v\\y 
— II'^IIei + ^ \^\ ■ 

□ 

Lemma 3.3. PyR^^m is a stationary solution to the Euler equations on Qr. 

Proof. Let am = ^VR'^m- Since ^pa^ and uj{am) are radially symmetric, so 
too must and uj{am) be. But then 

uj(am ■ Vam) = am ■ Va;(am) = 0, 

and thus am ■ Vct^ = Vp for some scalar field p. □ 



Remark 3.4. Equations ()3.4p through (j3.6p continue hold if the projection 
operator Pvr is replaced by the approximate projection operator \Jr, though 
a constant factor is introduced on the right-hand sides of Equation ()3.5p and 
Equation (j3.6p . Equations (j3.4p through (j3.6p also hold with Py^ replaced 
by U/j and i/^ replaced by L^. This gives control not only on the i?^-norm 
but individual control on the L^-norm. See Remarks 15.31 and 15.61 
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We will use the operator Fvr in establishing the deterministic expanding 
domain limit in Section [5] and in constructing statistical solutions to (NS) 
in Section [71 For solutions to the Euler equations, we will need the following 
approximate truncation operator of Definition 13. 51 (or we could use projection 
into Y{nn)). 

Definition 3.5. Define U/j : Y ^ ^(^r) by 

Uij(a„ + v) = (t^Iur + Trv, 
where : Yq Y{^}^) is the operator in Lemma 4.2 of [9j. 

Because is a disk, fTmlrj^ is in Hr and so also in Y(Qfi). This is why 
U/j: Y —f Y{Qji). It is also why if we define Prr be the restriction to 
followed by projection into Hr that for all u in E, 

PhrU = cf{u) + Phr{u - a{u)). 

It follows from Lemma 4.2 of [9j that for all u in E, 

\\PhrU - u\\ ^2(Q^) ^ as ii ^ oo. (3.8) 

Actually, Lemma 4.2 of [9] applies to an approximate projection operator 
into Hr, but we are using, as in in [lOlEI' the fact that projection into Hr 
gives the closest vector field in Hr to the vector field being projected. 

4. Weak deterministic solutions 



Definition 4.1 (Weak Navier-Stokes Solution). Given viscosity > 0, ini- 
tial velocity uq in Hr, and forcing / in L^^^([0, oo), //r), u in L?'{[^,T]]Vr) 
with dtu in L^([0, T]; V^) is a weak solution to the Navier-Stokes equations 
on ^Ir if n(0) = uq and 

(NS) / dtu-v+ (u- Vn) -v + u Vu-Vv = (uj) 
J^r J^r Jqr 

for almost all t in [0, T] and for all v in Vr. A weak solution on is defined 
for no in Em with u lying in L'^{[0, T];H'^) and dtu in L'^{[0, T];V'), and with 
(NS) holding for all v in V. 

Definition 4.2 (Weak Euler Solution). Given an initial velocity uq in Y{il.R) 
and forcing / in Lf^^i[0, (X)),Hr), u in L~([0, T];HRr\ H^nR)) with dtu in 
L^([0, T]; V^) is a weak solution to the Euler equations if u{0) = uq and 

(E) / dtu-v+ {u-Vu)-v = (/, v) 

J J ^R 

for almost all t in [0,T] and for all v in Hr n H^{Qr). A weak solution 
on is defined for uq in Y^ with u lying in L°°([0, T]; Y^) and dtU in 
L2([0, T]; V), and with {E) holding for all v in V. 
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Note that the test functions always have finite energy, even for solutions 
in the whole plane. (Test functions in Em would be too large to define the 
integrals involving the nonlinear terms in (NS) and (E).) 

Given a solution to {NS), there exists a distribution p (tempered, if the 
solution is in the whole plane) such that 

dtu + u-\/u + Wp = uAu + f, (4.1) 

equality holding in the sense of distributions. This follows from a result of 
Poincare and de Rham that any distribution that is a curl-free vector is the 
gradient of some scalar distribution. 

Given a solution to {E) , there exists a pressure p such that 

dtu + u-Vu + S/p = f, (4.2) 

but we can only interpret p as a distribution when working in the whole 
plane. Otherwise, we must view dtu + u ■ Vu as lying in H~^{ilii) and p as 
lying in L^(r2ij). (Equation (j4.2p follows, for instance, from Remark 1.1.9 p. 
14 of [TT].) 

In both Equation (j4.ip and Equation (j4.2p the pressure is unique up to the 
addition of a function of time. We resolve this ambiguity on by requiring 
that Jq p{t) = and on by requiring that p(t) lie in L?'{VLf{) for almost 
all t in fo,r]. 

In referring to solutions on 17 we will say solutions for R in [l,oo) and 
in referring to solutions on we will say solutions for R = oo. 

Theorem 4.3. (1) Assume that uq is in Em H . There exists a unique 
weak solution {u,p) to (NS) in the sense of Definition ^ . 1\ with initial velocity 
Uq for R = oo and initial velocity Py^uo for R in [1, oo), with 

u-amGL^i[0,T];HR), Vn G L~([0, T]; ^^(f^^)). (4.3) 

Moreover, there is a bound on each of these norms that is independent of R 
in [1,00] and that depends continuously on \\uo\\^i. For R = 00, if uq is in 
EmDH^ anduj{f) is in L'^{[0,T]; LPo n L°°) then uj{u) is m L~([0, T]; L^" n 
L~). 

(2) Assume that uq is in Y„ and that uj{f) is in L^([0,T];LPo n L°°). 
There exists a unique weak solution {u,p) to (E) in the sense of Defini- 
tion \4-^ with initial velocity uq for R = 00 and initial velocity U rUq for R 
in [1, 00). We have, 

u-ameL°°{[0,T];HR), V« G L-([0, r]|L2(0^)), 

u € L-([0, T] X Qn), u G C([0, T] x Q«), 

dtU G L-([0,r] ;//«), Vp G L-([0, T]; L2(17«)), 
a;(n) GL°°([0,T];LPonL°°(17fi)), 

and there is a hound on each of these norms that is independent of R in 
[1,cxd] and that depends continuously on \\uq\\y- 

Proof. These results are standard for R < 00, except for the independence 
of the norms on R. The independence of the first norm in Equation ()4.3p 
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follows for solutions to (NS) from the energy inequality in Equation ([52 
along with Equation (j3.4p : for the second norm it follows from adapting 
slightly the proof of this same fact for finite energy in [9]. The stronger 
bounds for solutions to {E) follow from the vorticity equation for (E). The 
results for R = oo are a minor modification of the same results for finite- 
energy: see, for instance, [2]. □ 

Definition 4.4 (Solution operators). Fix //j in L^([0, oo); i^/j) and write / 
for /oo- Let Sji{t) be the solution operator for (NS) on Or with S = Soo, 
and let Sii{t) be the solution operator for (E) on 17 with S = S^o- 

In application, we will often start with / in L^([0, oo); ff) or even time- 
independent f in H and let fji = Pjj^f. 

Sji{0) is the identity operator, as is Sr{0). For all t > 0, Sr maps 
Hr Vr and Vr Vr and S maps Em^ ErnHH^, E^ DH^ EmH H\ 
E ^ E\ and ^ For all t >0,Sr maps Y{nR) Y{nR) and S 
maps Ym — Ym and Y ^ Y. Each of these maps is continuous. 

Observe that SR{t)PvRUo = u{t) for R in [1, oo) and S{t)uQ = u(t) for 
i? = oo in part (1) of Theorem 14.31 while SR{t)'URUQ = u{t) for R in [1, oo) 
and S{t)uo = u{t) for i? = oo in part (2). 

5. Deterministic expanding domain limit 

First we establish the basic energy equality for deterministic solutions to 
(NS) in ^R and in all of M^. In all of M^, the energy is not finite, so we 
need to subtract am from the velocity to produce an "energy" equality. To 
make these estimates uniform over R in [l,oo], we need, then, to subtract 
am from the velocity for i? < oo as well. Actually, it will be slightly more 
convenient to subtract 

Wm = PyjjO-m for R G [1, oo), am = CTm for R = oo 

instead, but because Wm — > o"m in the H^[^R)-novm as i? ^ oo by Equa- 
tion (|3.4|) . this amounts to the same thing. 



Theorem 5.1. Let u be a solution to (NS) as in Definition \4.1\ with initial 
velocity uq in Hr for R < oo and uq in Em for R = oo. Then 



ft 

(n - CT„)(t)||^2(f^^) +2u I \\\7{u - am)\\l2^Q^) 





t 

JQr 



Wo - 0"m|lL2(f^^) - 2 / / {{u - am) -^CTm) ■ {u - am) (5.1) 



2u / Vam ■ V('u - am) + 2 j fR - {u- am) 
Jo JQ,R Jo Jn 



12 



JAMES P. KELLIHER 



and 

|2 , ,. / iiv7/„, ^ M|2 



\\{u-am){t)\\L2^n^^+l^ / l|V(u-CJ„)||^2(f^^) 

•^0 (5.2) 
< (iKo - + Cm'ut + ||/|lii([o,];L^)) e(^™+i)*, 

where C depends only on || Vfii |[^2nj;^oo • 

Proof. Assume first that R < oo. Using u — Wm, which is in Vr for all t > 0, 
as a test function in (NS) gives 

/ dtu-{u-am)+l (u-Vu) ■ (u-am) + 1^ Vu-V^u-Wm) 

f ■ (U-Wm)- 



R 



But, 



dtu- {u- Gm) = I dt{u- am) ■ (u- am) = l:^\\u- a. 



and 



{u ■ Vu) ■ [u - am) = {u- V(u - am)) ■ (u - am) 
+ (u- Warn) ■ (u - am) 

= -/ U ■ \/ \u - Wml"^ + {{u - Wm) ■ '^O'm) ■ {u - am) 

) ■ 

The first integral in the right-hand side above is formally zero because 
divii = and n • n = on d^R. More properly, we first observe that the 
integral is finite. This is because at time t in [0, T] both u and u — Wm are 
in H^{^r) and so in L^{Qr) by Sobolev embedding. But |V \u — Wm\ \ < 
2\u — Wm \ |V(u — Wm)\, and applying Holder's inequality gives the finiteness 
of the integral. Approximating by smooth functions and using the dominated 
convergence theorem shows that the integral is zero. (This is the approach of 
Lemmas II. 1.1 and II. 1.3 p. 108-109 of [H].) Since by Lemma [STSl Wm ■ ^Wm 
is a gradient, the last integral above vanishes. We conclude that 



(u • Vu) ■ {u - am) = {{U- CTm) " VcJm) • {u - am)- 

For the final term, we observe that 
/ Vu-V{u-Wm) 

JflR 

= V{u -Wm) -Viu -Wm) + VWm-V{u-Wm)- 
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Combining all these equalities gives 



2di 



\U - 0-m|li2(Q^) + I' ||V(U - 0-m)|li2(Q^) 



((n - ajn) ■ Vdm) • {u - am) 

Jr- {u- Wm)- 



+ 



Integrating in time gives Equation (j5.ip . Also, we can bound the right-hand 
side by 

llVa, 



'mil ijOo 

+ WfR 
< Cm \\u - 



I — Il2 

\u — a. 



'm|lL2(Q^) , , 2 

where we used Young's inequality and a 
d 
~dt 



V^iSlR) + ^ l|V'^m|lL2(R2) ||V(n - 0-m)||L2(n^) 
U - 0'm||2,2(Q^ 

1 . „ „9 ^ .,,2 



/filli2{Qfl) + ^ llV(n-cj^ 



m,o\. Thus, 



I — l|2 



+ zy||V(n-cJm)||i2(n^) 



< Cm^v + ||//?||i2(s^^) + C ||m - cj. 



^ ||2 



(5.3) 



Integrating in time and applying Gronwall's inequality gives Equation (j5.2p . 

The energy argument above works equally as well when R = oo with 
the exception of the term (1/2) f^2 u ■ V \u — Wmf , which must be handled 
slightly differently, because at time t in [0, T] we no longer have u in i?^(Ojj), 
only in i^^(R^). So we divide the integral in two, writing 

2 Jr2 



li • V lli — CTr; 



+ 



V7 I — |2 

V U - am i 



The first integral is finite and, in fact, zero, using the same reasoning as 
with the similar term for R < oo. The vector field a^n is in L°^(M?) 

andV|u — ami is in L^(M^) since |V |ii — am| | < 2 |u — am| |V(u — CTm)|; 
hence, the second integral is finite. We then have 

1 

2 

1 
2 



— V7 I — |2 
(Tm, • V n — (Tm, 







= lim 




R^oo 





(Jr. 



V7 I — |2 

V U - CJm, 



lim 

R^oo 



(o-r 



\ I — |2 
n - dm I 



□ 



This integral vanishes, since am • n = on OQr. 

Corollary 5.2. Let f lie in L^([0, 00); ff) and let Jr = PhrJ- Let uq be 

in Em n and let u be a solution to (NS) as in Definition \4- 1\ with initial 
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velocity Py^uo when R is in [l,oo) and initial velocity uq when R = oo. 
Then for a constant C independent of R and uq, 

\\{u-am){t)\\l2,Q^+l' I ||V(U - CTrn)|li2(Q^) 

•^0 (5.4) 

Proof. Apply Theorem 15.11 with initial velocity Pi/jjiio = ^Vr{uq — Om) + 
Pv^Cm and use Equation (j3.5|) . □ 

Remark 5.3. Were we to use an initial velocity of Uijuo instead of Pv^uq 
in Corollarv 15.21 we could replace H^io — '^mllyj; on the right-hand side of 
Equation (|5.4[) with C\\uo — UmWHR- See Remark 13.41 

We can control the decay of the tail of solutions to {NS) at time t based, 
ultimately, on the their decay at time zero: 

Lemma 5.4. For all uq in E^, 

||5(t)uo -o-(5'(t)no)||icx,([o,r];L2(Qg)) ^ as i? ^ oo 

and 

\\S{t)uQ - f^(5'(t)uo)||i2([o,r];jyi(f2g)) ^ as i? ^ oo. 

Proof. This is a minor adaptation of Lemma 7.1 of [9j to account for infinite- 
energy, and follows by a standard argument. □ 

Theorem 5.5. Assume that uq lies in ¥} and let uji(t) = Sji{t)Y'Yj^UQ and 
u{t) = S{t)uo. Then 

\Wr ~ '^\\l^{[o,t];H'^{Qii)) ~^ ^ as R ^ oo, (5.5) 

||V(ni?-u)||^2([o,T];L2(Cfl)) ^0 asR^oo, (5.6) 

and 

\\F{t,u) - i^i?(i,^^fl)|lL2([o,T];y^(nfl)) ^ as ^ oo. (5.7) 

In addition, the supremum over all uq in any bounded subset of and over 
all R in [1, oo] of each of the quantities, 

\\UR - f^(^)l!L°°{[0,T];L2{nfl,)) ' II ^""R II L2{[0,T];L2{nfl)) ' 
||i^fi(i,'"i?)|lL2([o,T];V^{Qfl)) 

is finite. 

Proof. The first two bounds in Equation (j5.8p follow from Equation (j5.4p . 
Equations (|5.5p and ()5.6p follow from Theorem 8.1 of [9] extended to infinite 
energy solutions using Equation (j5.4p in place of the standard finite-energy 
energy bounds. 
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We now prove Equation (j5.7p . We have, 

\\F{t,u) - FR(t,i^K)||y/(n«) < WMt) - 

+ \\Bu{t) - BRUR{t)\\y^(^^^^ + 11/ - • 

Let V be in Vr with H'^llv^ = 1- Then 

I {Au{t) - ARURit), v)\=u\ {Au{t) - AuR{t), v) \ 

= u\{Vu{t) - VuR{t),Vv)\ < V \Vu{t) - Vnfi(t)||^2(^^) \\v\\y^ . 

Thus, 

\\Au{t) - ^i?(ni?,)(t)||i2([o,T];y^(nB)) ^ ^ H'^^ " V'"RllL2([o,r];L2(Qfl)) > 

which vanishes as i? ^ oo by Equation (|5.6p . 
For the nonhnear term, 

\(Bu{t) - BRUR{i),v)\ = \{u{t) -Vuit) - UR{t) ■VuR{t),v)\ 

= |(div(n(t) (g) u{t) - UR{t) ® UR{t)), v)\ 

= \ {u{t) ® U{t) - UR{t) (g) UR{t), Vv)\ 

< \\u{t)0u{t) - nij(t) (g)nij(t)||^2(n^) \\v\\y^ . 

But, 

\\u{tyu{ty - UR{ty o URityw^^^^^^ 

< \\n{ty{u{ty - UR{ty)\\^,^^^^ + \\{n{ty - UR{ty)nR{ty\\^,^^^^ 

< c \\u{ty\^i^^^^ \\u{t) - UR{ty\^2^^^) 

< C Wum'/,]^^^ l|Vn(t)||i/'(^^) \\u{t) - URm'/,]^^^ 

X \\v{u{t)-uRm\'/,\^^^ 

where we used Ladyzhenskaya's inequahty (which gives no dependence on R 
for the constant C). Thus, 

\\Bu{t) - -Bij(n/j)(t)||^2([o,r];V^(nB)) 

< C \\'^u\\LmO,T];L^nR) 11^(^(0 - ^fl(0)llLi([0,T];L2{Qfl) 

Ct ||Vn||^2([o,T];L2(f7^) - UR{ty\\^2(^[o^T];LHnn) > 

which vanishes as i? ^ oo by Equation (|5.6p . 
For the forcing term, 

11/ - fR\\L^{[0,T];V^{nR)) - 11/ ~ /«llL2([o,T];//fl(f7fl)) ' 

which vanishes as i? ^ oo by Equation (|3.8p . 
From these bounds. Equation (j5.7p fohows. 

It remains to estabhsh the last bound in Equation (jS.Sp . But this follows 
from an argument similar to that we just made to prove Equation (|5.7p . 
using the first two bounds in Equation (|5.8p . □ 
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Remark 5.6. Together, the Hmits in Equations (jS.Sp and (j5.6p are called the 
expanding domain limit in [9] . Most of the estimates involved in establishing 
these limits require only that the initial velocity lie in H (or E for the infinite- 
energy extension). The key exception is that the regularity of the pressure 
is insufficient to complete the argument unless the initial velocity is in V (or 
for the infinite-energy extension). 

Had we used in place of Py^ in defining the initial velocity, the ex- 
panding domain limit would still hold (indeed, this is how the limit was 
established in [9] ) . An advantage of using is that the resulting bound on 
the rate of convergence is slightly improved, since i?^-norms are replaced by 
L^-norms in certain constants that appear in the bound. But this is unim- 
portant in our use of the limit, so we preferred to use Pvr, since it has a 
more natural definition. 

Corollary 5.7. For all uq in E^, 

\\SR{t)PviiUo - Pyii5'(t)uo||^oc([o,T];iTji) ^0 as oo (5.9) 

and 

\\SR{t)PvRUo - PyH'S'(*)'"o|lL2([o,T];Vii) ^ as ^ oo. (5.10) 
Proof. By Equations ()5.5p and ()5.6p 

^lirn^ \\SR{t)Pvj^Uo - S{t)uo\\^2(^[f)^T];Va) ^ as E ^ oo. 
But by Equation (i3l4]l . 

\\PvRS{t)uo - 5'(t)no||^2([o,T];yfl,) 

< C \\S{t)uo - a(5(t)no)||i2([o,T];i?i(Q^\n,/,)) + T'^' \m{u)\ /3(i?) 

< C \\S{t)uo - cT(5(t)no)||i2([o,T];HHng/,)) + T'^' \Hu)\m), 

which also vanishes as ii ^ oo by Lemma |5.4[ Equation (j5.10p then follows 
from the triangle inequality. The proof of Equation (j5.9p is similar. □ 

Remark 5.8. It is only in the proof of Corollarv 15.71 where we directly use 
the uniform decay over time of the tail of the velocity for solutions to {NS) . 
It was, however, already used in the extension of the expanding domain limit 
from finite to infinite energy energy alluded to in the proof of Theorem 15.51 

Lemma 5.9. For all u in E^. 

\\F{t,u) - FR{t,Pvjiu)\\y^f^^^^ as R ^ oo, 

and the supremum over any bounded subset o/E^ and over all R in [1, oo] of 

\\FRit,PvRu)\\y^^^^^ (5.11) 

is finite. 

Proof. The proof is the same as that of Equation (j5.7p . with no need to 
introduce the L^-norm over [0,T], and using Lemma 13.21 in place of the 
bounds in Equations (|5.5p and (|5.6p . □ 
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6. Definition of infinite-energy statistical solutions 

Following [6j p. 264-265, we define a statistical solution to (NS) on Qr, first 
defining the space of test functions. 

Definition 6.1. The space Tr of test functions, R < oo, is the set of all 
functions <I>: Hji — > M such that <I>(n) = (j){(u,gi), . . . , {u,gk)) for some (p in 
C^(R'') and some gi,. . . ,gk in Vr. The Frechet derivative of such a <I> is 
given by 

k 

^'(u) = ^djcl){{u,gi), {u,gk))gj, 
i=i 

which lies in Vr since each gj is in Vr. When i? = oo, we require that each 
oi gi, . . . , gk he compactly supported in M^, so that $: E — > M, and we also 
write T for T^q. This is the same class of test functions as for homogeneous 
solutions in the whole space (Definition 2.3 p. 278 of [6]). 

Observe that because each djcj) is bounded, 

\\<!>'{u)\\^^<C{<!>) (6.1) 

for all u in Vr. 

For t > and u in Hr let 

FR{t, u) = fR{t) - uAru - Br{u), 

where Ar is the Stokes operator and Br is the classical linear operator 
associated with the nonlinear term in (NS) on Qr. (See, for instance, p. 38 
of [6].) We also write F for F^o, A for Aoo, and B for B^q. 

For R = oo, we will assume for simplicity that / is in Lf^^{[0, oo); H); that 
is, we do not allow infinite forcing. 

Definition 6.2 (Statistical solution to (NS)). Assume that /hq is a Borel 
probability measure on Hr. Then a family, 

of Borel probability measures on Hr, R < oo, is a statistical solution to 
(NS) (SSNS) on $1^ if each of the following is satisfied: 
(1) For all $ in Tr and all t > 0, 

Hr JHr 



+ 11 {FR{s,u),^'{u))d^Xs{u)ds. 
Jo Jhr 



(2) For all t > 0, 



/ \\u\\\2 dutiu) + 2v / / II Vu 11^2 d/is(ti) ds 
J Hr Jo JHr 

/ {f{s),u)dns{u)ds+ \\u\\l2 d^oiu)- 
Jhr Jhr 
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(3) The map 



t 1-^ cl){u) djitiu) 
Jhr 

is measurable for all i > and all (/» in C^(Hji). 
(4) The map 




lies in L-([0,oo)). 



(5) The map 



t I— > / ||Vu||^2 dfifiu) 
Jhr 



lies in L;^^^([0, oo)). 



When R = oo, we make two changes in the definition. First, we replace 
by E throughout. Second, the energy equality in (2) is replaced by 



where a{u) is defined in Equation (j2.7p . 

The following is from Theorems 1.1 and 1.2 Chapter V of 
Theorem 6.3. Let fiQ be as in Definition W.i^ . R < oo, with kinetic energy 



and assume that f lies in L^^^([0, oo); There exists a SSNS, fi, as in 
Definition \6.2l If the support of fiQ is {Hn^Vn) -hounded as in Definition \2. 1\ 
(meaning that the containment in Equation l[6.3\) holds for t = 0) and f 
in Hn is time-independent then = Sii{t)^Q for all t > is a SSNS. 
Furthermore, this solution is the unique SSNS satisfying Equations i6.2\) 



(2') For all t > 0, 
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through 

1 1— > / ip{u) dfit{u) is continuous on [0, oo) for all ip in C{H'^), (6.2) 
Jhr 

supp^f Q Vr: \\u\\fj^ < m| for all t>0 for some M, (6.3) 
/ ^{t,u)d^it{u)= j ^(0,M)d^o(^t) 

J Hr J Hr 

+ f j [^',{s, u) + {F{u), u))] dfxsiu) ds. (6.4) 

Jo JHr 

H"^ is the space Hr in the weak topology. In Equation \6. 'J\) . equality holds 
for all Frechet-differentiahle continuous real-valued functions on [0, oo) x Vr 
(see the discussion following Equation V.1.16 in [6j for more details). 

For statistical solutions to (-E), we consider only solutions in the whole 
plane. For solutions to {E) there is no term involving the Stokes operator, 
so we define 

F{t,u) = f{t)-B{u). 

Definition 6.4 (Statistical solution to {E) in the plane). Assume that 
is a Borel probability measure on E. A statistical solution to the Euler 
equations (SSE) on E satisfies all the properties of a SSNS in Definition 16.21 
for R = oo except that the terms involving in property (2') are eliminated. 

7. Construction of Navier-Stokes solutions 

Let S{t) be the solution operator on E as in Definition I4.4[ Given that we 
expect the analog of Theorem 16.31 to hold for infinite-energy solutions in M^, 
we would expect that 

fit = S{t)fio (7.1) 

is the unique SSNS associated to the initial measure fiQ if we assume that 
the support of the initial Borel probability measure /io is (E,Ei)-bounded as 
in Definition 12.11 We show that this is, in fact, the case. Our approach will 
be to use the SSNS on Qr and take a limit as ii ^ oo in a careful way to 
demonstrate that fit is a SSNS on all of M^. 

We start by defining the initial probability measure fiQ on Hr by 

for all Borel measurable subsets E of Hr. Then fi^ is a probability mea- 
sure, for ij,(^{Hr) = iJ,o{Py^HR) = fio(E) = 1. Since we are treating initial 
probability distributions supported on E^, we use projection into Vr. When 
working with SSNSs as weak as those of Definition 16.21 projection into Hr 
would be used instead (though the limiting argument in that case is consid- 
erably more involved). 
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Similarly, we define the forcing term fn in by letting 

fR = PhJ, (7.2) 

where is projection into Hji. For simplicity, we assume that / is time- 
independent. Then 

ll/-/i?ll//H^Oasi?^oo (7.3) 

from Lemma 4.2 of [9j and the observation that projection into Hji gives the 
closest element in Hr. 

We let /x^ be the associated SSNS on so that, by Theorem 16. 3[ 

/if = SR{t)^,§, 

meaning that SR{t)fiQ-{E) = fiQ{S]^^{t)E) for any Borel measurable subset 
E of Hr. 

Let $ be in T as in Definition [HTTJ Because each gj is compactly supported, 
for all sufficiently large R, we can define a test function ^r in Tr by 

^r{v) ''^'(t>{{v,gi\nR),...,{v,9kK)) = '^{£Rv) (7.4) 

for all V in Hr, where £rv is extension by zero of v in Hr to all of M?. It 
follows that for all v in Hr, 

£r^'r{v) = ^'{£rv). (7.5) 

From now on, we always assume that R is sufficiently large that Equa- 
tion ([731) holds. 
For all n in E"*^, 

\^RiPvnn)-^iu)\ 

= \4>{(PvRU,9i\nii),-- ■ , {'PvRU,gkW)) - (l){{u,gi), {u,gk))\ 

< ||V(/)||^oo \ {PvrU - u, gi), . . . , {PvrU - u, gk)\ 

< UWc^ W^vrU - ^llL2(f7«) \ \W\\h + • • • + WgkWn ) 

Thus from Lemma [321 

^r{¥y^u) ^{u) as R ^ oo. (7.6) 
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Similarly, for all n in E^, 
|$^(Py^n)-cl>'(n)||^^ 



k 

< XI \dj4>{{'PvRU,gi), {'PvRU,gk)) - dj4){{u,gi),. . . , {u,gk))\ Ibjlly^ 

k 

^ C'^lJ'j{\{^VRU,gi),. . . ,(PvRU,gk) - {u,gi), . . . ,{u,gk)\) 
i=i 

k 

= C^fij {\(PvrU- u,gi), {PvrU - u,gk)\) 
i=i 

k 



.7 = 1 ^ 



<C>^fij [ ||Py^n-n||^ ( + - + \\9k\\H 



1/2 



J 

where is the modulus of continuity of dj(j). Thus by Lemma 

\\^'ji(PvRu) - <!>'{u)\\y^ ^ as ii ^ oo. (7.7) 

Before proceeding, we mention one logical simplification that we cannot 
make. It might seem reasonable to try to show that 

S{t)l^o = IhxL Pvr o Suit) o Pylno 

_i 1 (7.8) 

= hm /xo o o SR{t) o Pvr 

by showing that equality holds on any Borel measurable set E. We note, 
however, that if /xq is supported on a singleton set E = {uq} in E, uq nonzero, 
then Py^ o5'ij(t)^^ oPy^(S'^(t)'Uo) will in general never equal uq. Thus, the 
right-hand side above applied to E will evaluate to for all R, while the 
left-hand side will evaluate to 1. 

Observe that Equation (j7.8p is equivalent to saying that 



<^>{u) dfit{u) = lim / <i>{u) diy{^{u), (7.9) 
Jhr 

where ft = Py^^t for all test functions ^ (which are dense in the set of 
bounded continuous functions). We will prove instead that 

^(u) dfit{u) = lim / ^r{u) dfi^{u), (7.10) 
Je Jhr 

and that similar limits hold for the other integral in Property (1) of Defini- 
tion [6?2l thus circumventing our difficulty. In this weak sense, the expanding 
domain limit could be said to hold for statistical solutions. 

Theorem 16.31 continues to hold for i? = oo if we impose at the outset the 
condition that the initial velocity /io is supported in E^. This condition is 
required to allow us to take advantage of the expanding domain limit and 
related bounds from Section [5j The idea for proving existence is to first 
assume that the measure has bounded support in E^, apply the results of 
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Section [5l which require such support, then use the Unearity of properties 
(l)-(5) of Definition 16.21 to drop the boundedness assumption. This leads to 
Theorem EH 

Theorem 7.1. Let jiQ he as in Definition \6.2\ for R = oo, but supported in 
E^, and having "energy" 



I 

Je 



u\\^ dnQ{u) < oo. 



m- 



Assume that f lies in Lf^^{[0, oo); V). There exists a SSNS, fi, as in Defi' 
tion \5.^ If the support of fiQ is {K,E,^) -bounded as in Definition \2. 1\ ( meaning 
that Equation \7.12 ) holds for t = 0) and f in V is time-independent then 



fit = S{t)fio for all t > is a SSNS. Furthermore, this solution is the unique 
SSNS satisfying Equations { 7.11^ through Ii7.13\ ): 



t I— > / Lp{u) dfit{u) is continuous on [0, oo) for all ip in C{H^), (7.11) 
Je 

supp/ut C {n G ||n||i£ < M{t)} for all t > 0, (7.12) 



^'(t,u)(i/it(n) = / ^'(0,u)d//o(n) 
E Je 

t 



+ [ [ [f',{s,u) + {F{u),Kis,u))] dfis{u)ds. (7.13) 
Jo Jx 



In Equation (7.12), M is continuous on [0, oo). In Equation ( [7.i^| ), equal- 
ity holds for all Frechet-differentiable continuous real-valued functions on 
[0,oo) X V. 

Proof. Existence: Assume first that the initial Borel probability measure 
/^o is bounded in E^, meaning that 

supp/iQ C |n e E^ : < M} , for some M, 

and define fit by Equation (|7.ip for t > 0. By Equation (j4.3p it follows that 

supp^t C {-u G E^: IIuIIei < M(t)} (7.14) 

for some continuous function M. 

In Theorem l5.5( for initial velocity uq in E^ we defined UR{t) = 5j?(t)PvflUo 
and u{t) = S{t)uQ. In this proof we will be integrating over all initial veloc- 
ities in E and calling the initial velocity u, to agree with the notation of [6]. 
In this notation, Equation (jS.Sh and Equation (IS.lip become 

\\SR{t)'PvjiU - Cr(^i)|lLcx>([o,T];L2{Qfl)) ' ll^'S'i?(*)PVfl^llL2([o,T];L2(nB)) , 

||i^ij(t,SR(t)Pi/^n)|I^2([o,r];V^{Cfl)) > \\FR{t,Vvj,S{t)u)\\y^^^^^^ 

(7.15) 

are bounded on supp /io uniformly over all i? in [1, oo]. This will allow us to 
apply the dominated convergence theorem in several steps in our proof. 
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Using Equation (j7.4p . for all u in E^, 



$ ( lim SR{t)Pvuu) = $ ( lim SitSii{t)Pvuy) 
yR^oo J yR—^oo J 

= lim $ {eRSK[t)Vy^u) = lim {SR{t)PvnU) , 

/?— >oo R-^oo 

where we used the continuity of Thus we have, 

^{u)d^lt{u)= / ^{S{t)u)d^lo{u) 



(7.16) 



/ ^> ( lim S'/j(t)Py^ti ) duoiu) 
lim ^R{SR{t)PvRu) d^o{u) 
lim / ^R{SR{t)PvRu) dfio{u) 

JE 

lim / <^R{SR{t)v) dfi§{v) 
Jhr 

lim / ^R (n) dn^{u) 



giving Equation (|7.10|) . The first equality follows from Equation (|7.1|) . since 
the space of bounded continuous functions is dual to the space of Borel proba- 
bility measures. The limit in the second equality follows from Theorem 15.51 
The third equality follows from Equation (|7.16p . The fourth equality fol- 
lows by the dominated convergence theorem, since ^r is uniformly bounded 
over R in [l,oo] and fiQ is a finite measure. The fifth equality follows from 
Lemma [7.2[ The sixth and final equality follows in the same way as does the 
first. 

This shows that Equation ()7.10p holds for all i > 0, so if we can show that 

{F{s, u), ^'{u)) dfis(u) ds 

t (7-17) 

= hm / / {FR{s,u),<i>'ji{u))dfif{u)ds 
R^ooJo Jhr 

then we will have established the first property of Definition 16.21 for /x. 



Je 
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Toward this end, 

/ {F{s,u),^'(u))dfis{u)ds 
Je 



iF{s, S{s)u),<i>'iSis)u)) dfio{u) ds 
f I \\ul{FR{s,PvnS{s)u),^'{S{s)u))d^io{u)ds (7.18) 
lim {Fr{s, PvnS{s)u), ^'RiPvnS{s)u)) dfio{u) ds 



Je 



Je -R^oo 

= lim /* / {FR{s,Pv^S{s)u),<^'R{Pv^S{s)u))dfio{u)ds. 
R-^°°jQ Je 

The first equality follows from Lemma [7.3i The second equality follows from 
Lemma 15.91 and Equation (jT.lSp . since ^'{S{s)u) is bounded and compactly 
supported — and so also we can view the pairings as being in either the duality 
between V and V' or between Vr and V^. The third equality follows from 
Equation (j7.7p . The fourth equality follows from the dominated convergence 
theorem using Equation (j7.15p . 

We would like to commute the roles of the projection operator and the 
solution operator in the right-hand side of Equation (|7.18p to allow us to 
apply Lemma 17.21 To do this, we estimate, 

D{s,u) = \{Fr{s,Pv^S{s)u),<^>'r{Pv^S{s)u)) 

- (FRis, SRis)Pv^u), <l>'RiSRis)Pv^u))\ 

< \{Fr{s,Pv^S{s)u) - Fr{s,Sr{s)Pv^u),^'r{Pv^S{s)u))\ 

+ \{Fr{s, Sr{s)Pv,u), ^'r{Pv^S{s)u) - ^'R{SR{s)Pvnu))\ 

< \\FRis,PvnSis)u) - FRis,SRis)Pv,u)\\yJ^'R{PvnS{s)u)\\^^ 

+ \\Fr{s, SR{s)Pvnn)\\y^ \\^'R{PvnSis)u) - ^'r{Sr{s)Pv^u)\\^^ . 

Letting 

h{R,u,s) := \\<S>'r{Pv^S{s)u) - ^'R{SRis)Pv,u)\\^^ , (7.19) 

we have 



r / Dis, 

Jo Je 



< 



u) dfio{u) ds 

C f [ \\Fr{s,Pv^S{s)u) - FR{s,SR{s)Pv^u)\\y, dfio{u)ds (7.20) 

JO JE " 

+ / \\FR{s,SR{s)PvaU)\\y, h{R,u,s)dfio{u)ds. 

Jo JE ^ 
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Applying the Cauchy-Schwarz inequality, we can bound the last term by 

\\FRis,SRis)PvRu)\\L2^[0^t];VI,) ll^(^'^'-)llL2([0,t]) d^u) 



E 



< c 



||/i(i?,n,-)|lL2([o,i]) d^io{u), 



using Equation (j7.15p . 
By Equation (173]) . 

h{R,u,s) = \\^'iiiPvnS{s)u)-<S>'j,iSR{s)Pvnu)\\y^ 

= \\<^'{£RPvr,S{s)u) - <^'{SRSR{s)FvRu)\\y . 

By Equation (j6.ip . < Co for all v in Hr, for some Co independent 

of R. Thus, h{R,u,-) < 2Cq. Also, because ^' : H ^ V is continuous it 
follows from Equation (j5.9p that h{R,u, •) ^ as i? ^ oo for all u in E^. 
Hence, for all u in E^, 

ll^(-R>«r)lli2([o,t]) = ^ h{R,u,sfds 

by the dominated convergence theorem. But then also \\h{R,u, OIliaQot]) ^ 
2Cot^/^ and applying the dominated convergence theorem again gives 

WHR^u, OIlLado,*]) ^ as i? ^ oo. 

We conclude that the second term on the right-hand side of Equation (j7.20p 
vanishes as i? — > oo. 

For the first term in the right-hand side of Equation (j7.20p , 

\\FR{s,PvnSis)u) - FRis,SRis)PvnU)\\vi^ 

< \\Fr{s,PvrS{s)u) - F{s,S{s)u)\\y^ 

+ \\F{s,S{s)u) - FR{s,SR{s)PvRu)\\y, . 

Since dfj,Q(u) ds is a finite measure on [0,t] x E and the first term on the 
right-hand side is both bounded and vanishes as -R — > oo by Lemma 15.91 and 
Equation (|7.15p . after being integrated over [0,t] x E the first term vanishes 
as i? — > oo. The L'^{[0, t])-norm of the second term on the right-hand side is 
bounded on the support of hq by Equation (|7.15p and vanishes as — > oo 
by Equation (15. 7p : applying the Cauchy-Schwarz inequality followed by the 
dominated convergence theorem shows that 

/ / \\F{s,S{s)u) - FR{s,SR{s)PviiU)\\y, dfio{u)ds 
Jo Je ^ 

< t^^"^ jj\F{s,S{s)u) - FR{s,SR{s)Pvau)\\L2([o,t];Vl^) dMu) 

vanishes as i? — > oo. 
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We conclude that D(s, u) integrates to zero in the hmit as — > oo, mean- 
ing that 



giving Equation (j7.17p . completing the demonstration that property (1) of 
Definition 16.21 is satisfied for fj,. 

The other properties in Definition 16.21 follow more easily, using the dom- 
inated convergence theorem and the first two bounds in Equation (jS.Sp . 
Thus, we have established the existence of a SSNS for R = oo when the ini- 
tial probability measure has bounded support in E^. But we can drop this 
restriction by exploiting the inherent linearity in the definition of a SSNS, 
as done on p. 318 of [6]. This establishes the existence part of the theorem. 

Higher regularity: We now add the assumption that the support of 
is (E, E^)-bounded. Equation (j7.1ip and Equation (j7.12|) follow much as did 
properties (2) through (5) of Definition 16.21 Adding the assumption that / 
is time-independent. Equation (|7.13p follows for i? = oo in the same way it 
does for i? < oo. 

Uniqueness: The proof of uniqueness for i? < oo on p. 319-321 of [6] 
applies with the following two changes: First, in the Galerkin approximation 
we use a basis for E^ in place of the eigenfunctions of the Stokes operator 
(the spectrum no longer being discrete). Second, we use the energy bound 
in Equation (j5.2p for i? = oo in place of the bound involving the eigenvalue, 
Xm, of the Stokes operator. □ 

We used the following two elementary lemmas in the proof of Theorem 1 7. II 
Note that when we say that equality holds between two integrals when the 
integrands are only Borel measurable, we mean that either both integrals 
are defined and equal or that both integrals are undefined. We state the 
lemmas this way because in their application we do not always know a priori 
that the integrands are integrable. 

Lemma 7.2. For any Borel measurable function f on Hr, 



-K-^ooJq Jig 

Finally, using Lemma 17.21 and Lemma |7.3^ 






(7.21) 
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Proof. First observe that / o Pvr is Borel measurable on E because Pvr is 
Borel measurable (in fact, continuous) and / is Borel measurable, so the left- 
hand side of Equation (j7.2ip is well-defined. When / = xe, the characteristic 
function of a Borel measurable subset E of Hji, 



[ f(PvnU) df^oin) = M'PvlE) = 4iE) = [ f{v) dfi^ 
Je Jhr 



(v). 



Equation (|7.21|) then holds for simple functions by linearity, for nonnegative 
functions by the monotone convergence theorem, and hence for all Borel 
measurable functions. □ 

Lemma 7.3. For any function f that is Borel measurable on Hr, 
fiu)dti^{u)= [ f{SRit)u)dij.^iu). 

Hr J Hr 

When f is Borel measurable on X, 



f{u)dfit{u)= f{S{t)u)dfio{u). 
X Jx 

Proof. As in the proof of Lemma 17.21 equality holds for simple functions, 
then nonnegative functions, then all Borel measurable functions. □ 

8. Construction of Euler solutions 

We construct infinite-energy statistical solutions to the Euler equations by 
making a vanishing viscosity argument, using the infinite-energy statistical 
solutions to the Navier-Stokes equations that we constructed in Section [71 
For initial velocities as in Theorem 18.21 we have the following for SSNSs: 

Theorem 8.1. Assume that the support of the initial velocity fiQ for a SSNS 
with R = CO is bounded in Y as in Definition \2.1\ and that f is time- 
independent and lies in Yq. Then the SSNS also satisfies 

supp/it C {li e Y: IIuIIy < M(t)}, (8.1) 

for a continuous function M independent of u, and for all p in [po,oo], 

I \\uj{u)\\^p d^lt{u) < I \\uj{u)\\^p diXQ{u) + I ||a;(/(s))||^p ds. (8.2) 

JE JE Jo 

Proof. It is a standard result that 

MS{t)u)\\^, < Mu)\\^, + f ||^(/(t))||^, (8.3) 

JO 

for all u in Yq. To prove it for p = r/q in lowest terms, with r even, one 
takes the vorticity of Equation (|4.1|) . dtoo + u ■ Vto = uAto + uj{f), multiplies 
both sides by io^~^, and integrates over space and time formally to give 

Mml. +Pip-1) f \\u:^/^-^Vu\\l,^^,. = \\u:o\\%+p [\u;{f),u;P-'). 



JO 
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An approximation and smoothing argument is required to establish the 
equaUty rigorously, and it then follows for all p in [pojOo] by the conti- 
nuity of the norm as a function of p. Applying Holder's inequality gives 
Equation (j8.3p . 

Now assume that u = am + u is in Then dfU = dfV and Au = Av 
on nf, where Aam vanishes. Thus, the only additional complication in 
the argument above is the presence of the additional term (am ■ Vio,u)P~^) = 
(l/p)((Tm, Vi^P). But this vanishes formally by the divergence theorem, since 
am • n = on OQr, hence this term need not be accounted for in the ap- 
proximation and smoothness argument. 

Integrating Equation (j8.3p over E gives 

/ \\uj{S{t)u)\\^p dfio^u) < / \\uj{u)\\^p dfioiu) + / \\uj{f(t))\\^p. 
Je Je Jo 

(The last term has no dependence on u so the integral over E disappears, 

fiQ being a probability measure.) But |[lij(-)||j;^p : Y — > [0, cxo) is a bounded 

continuous function on supp/io so Equation ()8.2p follows from fit = S{t)fiQ, 

and Equation (18. ip follows from Equation ()8.2p . □ 

Theorem 8.2. Assume that ^uq is supported in Y with 

I^IIe dfio{u) < oo, 

E 

and assume that f is time-independent and lies in Yq. There exists a SSE, 
fj,, as in Definition \6.4\ One such solution is fit = S{t)iJ,o for all i > 0, where 
S{t) is the solution operator for the two-dimensional Euler equations in M? 
as in Definition \4.4\ Furthermore, if the support of hq is bounded in Y as in 
Definition \2.1\ and f is time-independent then this solution satisfies Equa- 
tion \8.1\) for some function M continuous on [0, oo) and Equation ^8.2\) . 

Proof. Assume first that the support of is bounded in Y. Define JL^ = 
S{t)fio, and let fi be the unique SSNS for R = oo given by Theorem 17.11 
with the same forcing and initial data as for the Euler equations. Let $ = 
(j){{u,gi), {u,gk)) lie in T. Then gi,...,gk are in V and 

k 

'^'(u) = '^dj(l){{u,gi),... ,{u,gk))gj G V, 
j=i 

k 

\/<P'{u) = Y,dj(t>{iu,gi),---,{u,gk))VgjeL^, 
i=i 

with 

\\^'{u)\\y<C, ||V^''(n)||^2 < C (8.4) 

for some constant C independent of u in E. 
Now, 

f ^{u)diit{u) = [ <i>{u)diio{u) [ [ {F{s,u),<^'{u))diis{u)ds 
Jy Jy Jo Jy 
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SO, using F = f — vAu — Bu and F = f — Bu, 




{F{s, u), $'(n)) dfig{u) ds 



JY 



$(u)d(/io - /^o)W + / / {F{s,u) - F{s,u),<i>'{u))dns{u)ds 

Jo JY 

{Bu, ^'{u)) d{fis - /"s)('") ds 




/O JY 

But fiQ =]Iq, so 



JY 

^{u)d{^iQ-JlQ){u) -V / / {Au,^'{u))d^is{u)ds 

Jo JY 



{Bu, ^'{u)) d{fis — iJ,s){u) ds. 




^{u) dipt - fJ't){u) = -f / / {Au,^'{u))dij.s{u)ds 

Jo Jy 



JY 



We have, 
and 



{Bu, <^'{u)) d{ns - 7lJ(n) ds. 
{Bu,^'{u)) = {u-Vu,^'{u)) 



{Au,<^'{u)) = -{Au,^'{u)) = (Vn,V^''(n)), 
since ^>'(n) is in V. Thus, 

(^{u) d{fjLt -'p-t){u) = -V I I {\/u,V^'{u)) dfis{u) ds 
rt 



JY 



We have. 



- [ [ {u-Vu,^'{u))d{fj.s -'Ps){u)ds. 
Jo Jy 

{Vu,V^'{u))dns{u) <C ||Vii||j^2 dfisiu) < C, 

Jy 



where we used Equation (|8.4|) followed by Equation (|8.2|) and the bounded- 
ness of fiQ in Y. The same bound holds when integrating against Jig. Thus, 



$(n) d{fit - ^^t){u) 

■t 



< Rvt + 



JY 



[u ■ Vu, ^'{u)) d{fis — ^J's){'^) ds 



(8.5) 



where R is proportional to the right-hand side of Equation (|8.2p , which we 
note increases with time unless there is zero forcing. 
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For any Borel measurable function G on 

G{u)d{fis —'Ps){'^) = / G{u)d^s{u) — I G{u)d'i2g{u) 
Jy Jy 

G{S{s)u) dfio{u) — / G{S{s)u) dpQ{u) 
Jy 

{GiSis)uo)-G(S{s)uo))dMuo) 

(GKs))-G(n(s)))d/io(no). 

In the last integral, we are defining u{t) and u{t) to be S{t)uQ and S{t)uo, 
respectively. These are the solutions to (NS) and (E) given the initial 
velocity uq. (The support of fiQ lying in Y insures that S{t)uQ is well-defined 
and continuous for /LiQ-almost all uq.) 
Thus, 

{u-Vu,^'{u))d{i^s - Tis){u) 

[{u{s) ■ Vu{s),^'{u{s)) - (u(s) • Vn(s),$'(u(s)))] d^o- 

Letting w = u — u, we have 
(n • Vn, <^'iu)) - {u ■ Vn, <^'(u)) = {u ■ Vw, ^'(u)) 
+ {u ■ Vu, - ^>'(m)) + (w ■ Vu, ^'{u)) 
= -(n • V^'{u),w) + {u ■ Vu, - ^'(u)) + {w ■ Vu, ^'{u)), 

so 

|(n(s) • V'u(s),$'(u(s)) - (n(s) • Vu(s), $'(u(s)))| 

< ||n(s)||ioo ||V«>'(u(s))||^2 

+ \\u{s)\\l^ \\Vuis)\\^, ll^'Ks)) - ^'(uis))\\^ 

+ \\w{s)\\H\\Vu{s)\\L,\\'^'iu{s))\\^^. 

Now, Equation (j8.3p holds for solutions to (E): it can be derived as for 
(NS) or by viewing (E) as a non-homogeneous transport equation for the 
vorticity. Since supp /xq is bounded in Y, it follows from Equation ()8.3p that 
u and u are bounded in the L°°([0,r] x M^)-norm uniformly over supp/zo; 
as is Vu in the L°°([0, T]; L^)-norm. This is discussed more fully in [2j or 
[8], where it is shown, moreover, that there exists a continuous function 
p : [0, oo) X [0, oo) — > [0, oo), nondecreasing in t, with p{0, t) = for all t > 0, 
such that for alH > 0, 

\\w{t)\\H<p{u,t). (8.6) 



(For sufficiently smaU i^t, p{v,t) = (Ci^t)(i/2)e --^^ 



-ct 
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Also, 

|$'KS))-$'(II(S))||^ 



k 



^ X] • • • , {u{s),gk)) - dj(t){{u{s),gi), {u{s),gk))\ \\gj 

Now, 

\{u{s),gj) - {u{s),gj)\ < \\w{s)\\j^ \\gj\\^ < pii^, s) \\gj\\^ , 
so since each djcj) is continuous, it follows that 

||<I>'(n(s)) - ^'{u{s))\\^ ^ as ^ uniformly over [0,r]. 
Combining all these facts shows that 



/ {u- Vu, d{ns - fJ's){u) ^ as ^ 

Jo Jy 



and hence that 



lirn I ^{u) dfit{u) 



Jo Jy 



(8.7) 



On the other hand, 



lim / ^{u) dfitiu) = lim / ^{S{t)u) dfio{u) 



lim <^{S{t)u) d^oiu) = / <^{S{t)u) duoiu) (8.8) 
u^o Jy 

^{u)d4lt{u). 

In the second equality we used the dominated convergence theorem. For the 
third equality, we used 

\^{S{t)u) -<^(S{t)u)\ 

= \(t>{{S{t)u, gi),..., {S{t)u, gk)) - ct){{S{t)u, gi),..., {S{t)u, gk))\ 

< II V^ll^oo \{{S{t)u, 5i), . . . , {S{t)u, gk)) - {(S{t)u, gi), (S{t)u, gk))\ 

< C\{{S{t)u-S{t)u,gi), {S{t)u-S{t)u,gk))\ 

< C\\S{t)u-S{t)u\\jj < Cp{v,t) ^ as ^ 0, 

the last inequality just being another way of writing Equation (j8.6p . Hence, 
the right-hand sides of Equations (jS.Th and (|8.8p are equal, establishing the 
first property in Definition 16.41 

Equations (jS.ip and ()8.2p follow as in the proof of Theorem 18.11 
As in the proof of Theorem 17. H we can drop the restriction that the 
support of jJiQ is bounded in Y by exploiting the inherent linearity in the 
definition of a SSE, as done on p. 318 of [6]. The remaining properties in 
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Definition 16.41 follow using the dominated convergence theorem in a manner 



Since the space T of test functions is dense in the space of all bounded 
continuous functions on E, it follows that /i — > /J as measures as ^ 0; 
that is, the vanishing viscosity limit holds for statistical solutions to the 
Navier-Stokes and Euler equations. 



The author was supported in part by NSF grant DMS-0705586 during the 
period of this work. 
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similar to what we did above. 



□ 



The proof of Theorem 18.21 shows that 
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